Symmetry energy softening in nuclear matter with non-nucleonic 

constituents 
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We study the trend of the nuclear symmetry energy in relativistic mean-field models with appear- 
ance of the hyperon and quark degrees of freedom at high densities. On the pure hadron level, we 
focus on the role of A hyperons in influencing the symmetry energy both at given fractions and at 
charge and chemical equilibriums. The softening of the nuclear symmetry energy is observed with 
the inclusion of the A hyperons that suppresses the nucleon fraction. In the phase with the admix- 
ture of quarks and hadrons, the equation of state is established on the Gibbs conditions. With the 
increase of the quark volume fraction in denser and denser matter, the apparent nuclear symmetry 
energy decreases till to disappear. This softening would have associations with the observations 
which need detailed discriminations in dense matter with the admixture of new degrees of freedom 
created by heavy-ion collisions. 

PACS numbers: 21.65. Ef, 21.60.Jz, 21.65. Qr, 13.75.Ev 



I. INTRODUCTION 

The nuclear symmetry energy of isospin asym- 
metric nuclear matter is not only important 
for understanding the structure of neutron- or 
proton-rich nuclei and the reaction dynamics of 
heavy- ion collisions, see, e.g., Ref. [1-3], but 
also plays a crucial role in a number of impor- 
tant issues in astrophysics, see, e.g., Refs. [4- 
6]. Recently, appreciable progresses have been 
achieved on constraining the symmetry energy 
at saturation and subsaturation densities either 
through the extraction based on astrophysical 
observations or in terms of terrestrial data [7- 
12]. However, the density dependence of the 
symmetry energy is still poorly known at supra- 
normal densities [3, 13-15]. Theoretical models 
predict diverse density dependencies of the sym- 
metry energy at high densities. Noticeably, quite 
different density-dependent trends of the sym- 
metry energy can be extracted from analyzing 
the FOPI/GSI data on the 7r~/7r + ratio in rel- 
ativistic heavy-ion collisions with various trans- 
port models [13-15]. In spite of this inconsis- 
tency, the theoretical uncertainty of high-density 
symmetry energy is regarded to be associated 
with the tensor force that originates from the ex- 
change terms [16, 17]. In the ladder approxima- 
tion, the exchange terms can be well treated in 
the Brueckner theory either in the relativistic or 
non-relativistic frameworks [18, 19]. In deed, the 
vacuum polarizations given by ring diagrams arc 
absent in the Brueckner theory. The inclusion of 
the ring diagrams is however very complicated. 
In this work, we do not carry on the tensor force 
that appears beyond the Hartree approximation 
but consider the non-nucleonic degrees of free- 
dom in the Hartree approximation. 



The new degrees of freedom considered here 
are hyperons and quarks that may appear in 
dense matter roughly around the density 2-4po, 
depending on the paramctrization of models [20- 
24]. Nuclear matter at this density domain can 
be produced via heavy-ion collisions, and it usu- 
ally includes the admixture of non-nucleonic de- 
grees of freedom. In the past, the symmetry 
energy effects on these phase transitions have 
been found to be very significant [23, 25, 26]. 
However, the effects of new constituents on the 
symmetry energy are seldom investigated. It is 
not the aim of this work to resolve the uncer- 
tainty of the high-density symmetry energy but 
to reveal the variation of the symmetry energy 
in phases mixed with these new constituents. 
Once the phase transition occurs, the system 
goes to the mixed phase that can be theoreti- 
cally constructed by virtue of the phase equilib- 
rium conditions, namely, the Gibbs conditions 
in this work. In the mixed phase, we define the 
symmetry energy according to the general ex- 
pression of the energy density that is different 
from that in pure nuclear matter. The paper is 
organized in the following. In Sec. II, we present 
necessary formulas for the nuclear symmetry en- 
ergy in pure hadron and mixed phases with the 
relativistic mean-field (RMF) framework. The 
construction of the mixed phase of quarks and 
baryons arc presented briefly. In Sec. Ill, the 
numerical results and discussion are given. At 
last, we give a brief summary. 



II. FORMALISM 

In the parabolic approximation, the energy 
per nucleon in isospin asymmetric nuclear mat- 
ter can be written as 
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where eo(/o) is the energy per nucleon in symmet- 
ric nuclear matter, the E sym (p) is the symmetry 
energy, and 5 = (p n — p p )/p is the isospin asym- 
metry. In RMF models, the energy density can 
generally be written as [27] 
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where = g*/m* , C p = g*/m*, C a = m* a /g*, 

E* = ^k 2 + m* N 2 with m* N = M* - g*a the 
effective mass of nucleon, kp is the Fermi mo- 
mentum, and Enon is the nonlinear meson self- 
interacting energy density specifically for the 
nonlinear RMF models [28-30]. Here, the pa- 
rameters with asterisks denote the density de- 
pendence which is given in specific models. For 
instance, in nonlinear RMF models, the cou- 
pling constants are density- independent, while 
the meson masses are density-dependent due 
to the nonlinear self-interacting terms. With 
Eq.(2), the symmetry energy in the RMF models 
can be derived as 
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For hyperon degrees of freedom, we consider 
only the A hyperon for simplicity. In this case, 
Eq.(l) still holds for hyperonized matter. The 
nuclear symmetry energy now reads 
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where kp is the nucleon Fermi momentum, pjq 
is the number density of nucleons, and ps = 
Pn + PA- This formula applies to the case of 
the given ratio pa/pb- The symmetry energy is 
now suppressed due to the factor pn/pb- On 
the other hand, as one source term of meson 
fields, the A hyperon has led to a moderate de- 
crease to the nucleon effective mass and Ep in 
the kinetic term. Together with the suppressed 
nucleon Fermi momentum in Ep, the suppres- 
sion of the kinetic term can be partially com- 
pensated. Nevertheless, the symmetry energy 
eventually turns out to be suppressed in either 
the baryon-density or nucleon-density profile. In 
chemically equilibrated and charge neutral mat- 
ter, the particle fractions are obtained from solv- 
ing coupled equations. In this case, the nuclear 
symmetry energy is given as 



E S ym — TjiPn ~ Pp) — i (5) 
10 pB 



where p n and p p are the neutron and proton 
chemical potentials, respectively. As the isospin 
asymmetry parameter S approaches to vanish- 
ing, Eq.(5) reduces to Eq.(4). 

After the hadron-quark phase transition oc- 
curs, hadrons and quarks coexist in a mixed 
phase. The construction of the mixed phase is 
based on the mechanical and chemical equilib- 
riums, namely, the Gibbs conditions which are 
given as [21] 
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p y , p u = p n /3 - 2p e /3, 
Ps = Pnj'i + Pe/3. (6) 



The pressures of nuclear and quark matter read 
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where B is the bag constant of the MIT 
model [31] and So is the rearrangement term 
induced by the density dependence of model pa- 
rameters [27]. For convenient narration, we do 
not include the A hyperons whose addition can 
be referred to Rcf. [23]. In actual calculations, 
we include the strange mcson-hyperon interac- 
tions. In terms of the quark phase proportion 
Y, the total density can be expressed as 

PB = jPQ + (l-Y) PH , (9) 

where pn is the baryon density on the hadronic 
level and pq is the quark density. Using Gibbs 
conditions, one can obtain the quark phase pro- 
portion Y. The total energy density and isospin 
asymmetry parameters arc written as 

£ = (1 - Y)£ H + Y£ Q , a = (1 - Y)8 H + YS Q , 

(10) 

with S H = {Pn - Pp)/pn and Sq = (p u - 
Pd)/{Pu + Pd)- The energy density in the mixed 
phase thus depends on the Y. In the parabolic 
approximation, the energy density can be ex- 
pressed as 

£/p B = e (pB,Y)+E? ym (p B ,Y)6 2 H 

+E% m (pB,YW Q . (11) 

Because the volume fraction depends on the 
isospin asymmetry, we limit the derivation of 
the symmetry energy in symmetric matter, 
namely a = 0. In this way, the nuclear symme- 
try energy is defined as Ef ym = ±d 2 (£/p B )/d5 2 H 
at S h = 0, and the definition of the quark 
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symmetry energy is similarly given as E^ ym = 
\d 2 {E/p B )/d8 2 Q at S Q = 0. Bridged by Gibbs 
conditions, the quark volume fraction is model 
dependent and relies on the MIT bag constant. 
As a result, similar dependencies on the model 
and bag constant can be delivered to the sym- 
metry energy. As seen from Eqs.(10) and (11), 
the nuclear symmetry energy may disappear as 
the quark fraction grows to be unity at high den- 
sities. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

We first display numerically the role of A hy- 
perons in affecting the nuclear symmetry energy. 
In Fig. 1, it shows the density profile of the sym- 
metry energy for various A fractions with RMF 
models SLC and SLCd [32]. The symmetry en- 
ergy for various A fractions is calculated in sym- 
metric matter at S = 0. It is shown in Fig. 1 that 
the symmetry energy is softened clearly with the 
increase of the A fraction. For chemically equi- 
librated and charge neutral matter, we see that 
the symmetry energy starts to soften once the A 
hyperons appear. In this case, the symmetry en- 
ergy at lower densities with the SLC is identical 
with that obtained with /a = in symmetric 
matter, while a small difference appears in re- 
sults with the SLCd, which is associated with 
the deviation in calculating the symmetry en- 
ergy at the clearly larger isospin asymmetry in 
SLCd. The usual case denoted in Fig. 1 means 
that the hyperon-nucleon interaction has simi- 
lar in-medium properties to that of the nucleon- 
nuclcon interaction (for details, see Ref. [24]). 
Shown in Fig. 2 is the symmetry energy for var- 
ious A hyperon fractions with the separable case 
whose density dependence of the hyperon poten- 
tial is different from that of the nucleon poten- 
tial [24]. Except for the chemically equilibrated 
and charge neutral case, results shown in Fig. 2 
and 1 arc almost identical. This can be elabo- 
rated by Eq.(4) because the nuclear symmetry 
energy in hyperonized matter is dominantly af- 
fected by the hyperon fraction. In the separable 
case, the hyperon fraction saturates at the cer- 
tain high density, and the hyperons disappear at 
very high density. Thus, in the density profile of 
the symmetry energy, a concave segment forms. 

With the increase of density, the hadron-quark 
phase transition may occur. In this work, quark 
matter, regarded as the free fermion gas with- 
out interactions, is described with the MIT bag 
model [31]. For a simple example, we first 
deal with the phase transition without hyperons. 
In a next step, we then include the hyperons. 
The mixed phase consists of high-density quark 
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FIG. 1: The symmetry energy as a function of den- 
sity for various A fractions in the usual case (For 
details, see text). The upper panel is for the re- 
sults with the SLC, while the lower panel displays 
the results with the SLCd. The label /3 denotes the 
chemically equilibrated and charge neutral matter. 
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FIG. 2: The same as shown in Fig. 1 but with dif- 
ferent density dependencies for the hyperon and nu- 
cleon potentials, denoted as the separable case. 



matter and low-density nuclear matter with the 
quark phase proportion Y being obtained ac- 
cording to Gibbs conditions. Here, we do not 
iterate the solution details which can be found 
in the literature [21]. The critical density of the 
hadron-quark transition is model-dependent and 
also depends on the bag parameter clearly. In 
Fig. 3, we depict the critical density function 
of isospin asymmetry for various bag parameters 
and RMF models. We can observe a few char- 
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FIG. 3: The critical density for the hadron-quark 
phase transition as a function of isospin asymme- 
try with various RMF models and bag constants as 
marked. 



acteristics in Fig. 3. First, the critical density 
obtained with the TM1 [29] is much more sensi- 
tive to the isospin asymmetry than other mod- 
els. The difference may be associated with the 
softening of the vector repulsion in TM1, while 
other models shown in Fig. 3 do not possess this 
softening. Second, the critical density increases 
clearly with the rise of the bag constant. Third, 
the softening of the symmetry energy just has 
a moderate effect on the critical density. The 
RMF model SLCd has a softer symmetry energy 
than the SLC. The softening of the symmetry en- 
ergy also occurs in the RMF model NL3w3 [30] 
as compared with the model NL3 [28]. One can 
see that the sensitivity of the critical density to 
the isospin asymmetry reduces generally due to 
the softening of the symmetry energy. The ex- 
ception occurs for models SLC and SLCd with a 
smaller bag parameter, see left panel in Fig. 3. 

Besides the critical density, the quark phase 
proportion Y also depends on the isospin asym- 
metry. In this way, the symmetry energy in the 
mixed phase obtained in symmetric matter can 
not simply be used to predict the properties of 
asymmetric matter because the quark volume 
fraction changes with the isospin asymmetry in 
asymmetric matter. Nevertheless, the symmetry 
energy obtained in symmetric matter is instruc- 
tive to exhibit its variation trend in the mixed 
phase. Shown in Fig. 4 is the nuclear symmetry 
energy as a function of baryon density with the 
bag constant B = (18QMeV) 4 . Apparent de- 
crease of the symmetry energy can be observed 
after the hadron-quark phase transition occurs. 
With the increase of density, the volume fraction 
of nucleons decreases, which causes a straightfor- 



FIG. 4: The nuclear symmetry energy as a func- 
tion of density in symmetric matter including the 
hadron-quark phase transition at high densities. 
The reflection point from rising to dropping corre- 
sponds to the critical density for each model. Above 
the critical density, plotted is the symmetry energy 
F H 

^ sym • 

ward reduction of the nuclear symmetry energy. 
As the nucleon volume fraction reduces to zero, 
the nuclear symmetry energy vanishes. 

For a smaller bag constant, the critical density 
is also smaller, and the decrease of the symmetry 
energy starts at lower densities for various mod- 
els. Shown in Fig. 5 is the symmetry energy with 
B = (lQOAleV) 4 . In this case, the critical den- 
sity is around 1.6-2.2po and the symmetry en- 
ergy vanishes at densities below 3po for all mod- 
els used in the calculation. For a larger bag con- 
stant, quite large difference in the critical density 
and vanishing density of the symmetry energy 
exists for different models. This can be observed 
by comparing Figs. 4 and 5. For a much larger 
bag constant, for instance, B = (200MeV) 4 , the 
divergence becomes more appreciable. Never- 
theless, once the hadron-quark phase transition 
occurs, the decrease of the symmetry energy is 
definite. 

For the quark phase, the quark symmetry en- 
ergy Ef which reflects the cost deviating from 
the flavor symmetric matter starts to have value 
above the critical density and increases with the 
rise of the quark volume fraction in the mixed 
phase, as shown in Fig. 6. After the quark vol- 
ume fraction quickly develops to be unity that 
is a value for pure quark matter, the quark sym- 
metry energy grows quite slowly with the den- 
sity, since without interactions only the kinetic 
energy contributes to the symmetry energy. 

With the inclusion of hyperons, the re- 
construction of the chemical equilibrium with 
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FIG. 5: The same as shown in Fig. 4 but with B = 
{imMeV) 4 . 
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FIG. 6: The quark symmetry energy as a function 
of density for various bag constants and RMF mod- 
els. The difference in the symmetry energy exists 
in the mixed phase for various RMF models, and it 
disappears in pure quark matter at sufficiently high 
densities. 



quarks results in a larger critical density. Shown 
in Fig. 7 is the nuclear symmetry energy in the 
hadronic and mixed phases for various fractions 
of A hyperons with RMF models SLC and SLCd. 
Here, the bag constant is B = (180Afey) 4 . The 
inclusion of A hyperons suppresses the symme- 
try energy in the hadronic phase, consistent with 
those shown in Figs. 1 and 2. In the mixed 
phase, the influence of A hyperons is not promi- 
nent due to the suppression factor (1 — Y) that 
decreases quickly. For other bag constants and 



FIG. 7: (Color online) The symmetry energy as a 
function of density for different A-hyperon fractions 
in hyeronized matter with the hadron-quark phase 
transition. The RMF models SLC and SLCd are 
adopted here and the bag constant is (l&OMeV) 4 . 



other RMF models, the conclusion is qualita- 
tively similar. Namely, the inclusion of A hyper- 
ons suppresses the symmetry energy in hadronic 
phase, and the decrease of the symmetry energy 
starts at a larger critical density. To save space, 
these numerical results are thus not displayed. 

We note that the softening of the symmetry 
energy in the mixed phase is mostly apparent 
because once the quark volume fraction can be 
identified at given densities the nuclear sym- 
metry energy would be extracted appropriately 
by singling out the effect of suppression factor 
(1 —Y). However, the determination of the Y 
is strongly model-dependent and far from ex- 
perimental feasibility. On the other hand, the 
decreasing factor (1 — Y) in the mixed phase 
largely suppresses the growth of the nucleonic 
density with the rise of the total density. If 
the hyperons are taken into account, the sit- 
uation becomes more complicated. Thus, the 
extraction of the high-density symmetry energy 
for pure nucleonic matter is not well grounded 
once the hadron-quark phase transition occurs. 
Most likely, the high-density symmetry energy 
extracted from the heavy-ion collisions would be 
as soft as that presented in this work as long 
as a detailed discrimination or calibration is not 
ready for dynamically evolutional matter. 

At last, it is worth mentioning that the ap- 
pearance of new degrees of freedom that is en- 
ergetically favored usually softens the equation 
of state, resulting in a significant decrease of 
the maximum mass of neutron stars. Recently, 
the pulsar J1614-2230 was identified rather accu- 
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rately through the Shapiro delay to have a mass 
2Mq [33] , which sets up a lower limit of the max- 
imum mass of neutron stars. This states that the 
nuclear equation of state should not be softened 
significantly even with the appearance of new 
degrees of freedom. A way out of this dilemma 
is to consider new forms of interactions for new 
degrees of freedom [24, 34-38] . For quark mat- 
ter, the imposition of interactions can stiffen the 
equation of state and hence increase the max- 
imum mass of neutron stars [36-38]. It would 
be interesting to investigate whether the interac- 
tions of quarks have an effect on the nuclear sym- 
metry energy. This deserves subsequent work 
and is however beyond the scope of the present 
work. 



quark degrees of freedom, we have constructed 
the isospin symmetric mixed phase according to 
Gibbs conditions using the RMF models and 
MIT bag model. The nuclear symmetry en- 
ergy obtained in the mixed phase reduces quickly 
with the rise of quark volume fraction. We have 
recognized that the specific softening depends 
on the parametrizations of models. Especially, 
it has a clear dependence on the bag constant of 
the MIT bag model. Nevertheless, we conclude 
that the effect of phase transitions is important 
on the symmetry energy, and for the experimen- 
tal extraction of the symmetry energy at high 
densities it is significant and necessary to take 
into account the effect of phase transitions. 



IV. SUMMARY 

We have studied the effect of A hyperons and 
quarks on the nuclear symmetry energy at high 
densities with relativistic models. The softening 
of the nuclear symmetry energy is observed ei- 
ther in chemically equilibrated matter or matter 
with an given A fraction. With the inclusion of 
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